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We find concrete evidence for a recentiy discovered form of intermittency, referred to as in-out 
intermittency, in botli PDE and ODE models of mean field dynamos. This type of intermittency 
. (introduced in occurs in systems with invariant submanifolds and, as opposed to on-off inter- 

' mittency which can also occur in skew product systems, it requires an absence of skew product 

, structure. By this we mean that the dynamics on the attractor intermittent to the invariant man- 

ifold cannot be expressed simply as the dynamics on the invariant subspace forcing the transverse 
' dynamics; the transverse dynamics will alter that tangential to the invariant subspace when one is 

far enough away from the invariant manifold. 

Since general systems with invariant submanifolds are not likely to have skew product structure, 
I this type of behaviour may be of physical relevance in a variety of dynamical settings. 

The models employed here to demonstrate in-out intermittency are axisymmetric mean-field 
dynamo models which are often used to study the observed large scale magnetic variability in the 
Sun and solar-type stars. The occurrence of this type of intermittency in such models may be of 
interest in understanding some aspects of such variabilities. 

u 

■ Dynamical systems that possess symmetries (and hence invariant submanifolds embedded in their 
state spaces) are of interest in a variety of settings. In many simplified models such dynamical systems 
have skew product structure. For an ODE model, if {x, y) parameterizes a phase space with an invariant 
manifold y = 0, we say the system has skew product structure ii x = f{x) and y — g{x,y), namely if 
the dynamics of x is independent of y. A great deal of effort has gone into the study of such skew 
J>. product systems with invariant manifolds, and these have thrown up a number of new and interesting 
fS| phenomena. 

CO In general, however, one would expect dynamical systems not to have skew product structure unless 

extra structure is present (for example if the transverse dynamics is always forced by the tangential 
Jli^ dynamics). In the absence of such extra structure it is therefore interesting to see what new types of 
dynamics can appear in systems with invariant submanifolds. One such novel type of dynamical be- 
haviour, in out intermittency, is discussed and analysed in detail in ||^ using a simple two-dimensional 
mapping. An important feature of this type of intermittency is that, as opposed to on off intermit- 
tency, it requires the absence of a skew product structure. 

In this paper we find concrete evidence for the occurrence of in out intermittency in both PDE and 
ODE models both in terms of phase space and also statistically. The models considered are examples 
K>" of axisymmetric mean field dynamo models which are often used in order to study the observed large 
KJi ^ scale magnetic variability in the Sun and solar-type stars. In addition to providing examples of in— out 
intermittency in PDE models, the occurrence of this type of intermittency in such models may be of 
d interest in understanding some aspects of solar and stellar variabilities. 
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I. INTRODUCTION 



Many systems of physical interest possess symmetries which in turn induce invariant submanifolds in their state 
spaces. A great deal of effort has gone into the study of the dynamics and intermittent behaviour of such systems 
near their invariant submanifolds (see e.g. [Q). A class of dynamical systems with invariant submanifolds have been 
shown to be capable of producing a number of novel modes of behaviour, including on-off intermittency, which occurs 
as the result of an instability of an attractor in an invariant submanifold . It manifests itself as an attractor whose 
trajectories get arbitrarily close to an attractor for the system in the invariant submanifold while intermittently 
making large deviations away. It can be modelled by a biased random walk of the logarithmic distance from the 
invariant submanifold 

Since the linearised behaviour near an invariant submanifold has a natural skew product structure (i.e. the linearised 
dynamics transverse to the invariant submanifold is forced by the dynamics within the submanifold) many such studies 
have tended to concentrate on systems that are of skew product type for simplicity, although it should be stated that 
on-off intermittency can be found in systems that do not have skew product structure. 

Moreover, bifurcation problems in such settings have tended to concentrate on normal parameters, i.e. parameters 
that vary the global dynamics without changing the dynamics within an invariant submanifold. In general, dynamical 
systems are not skew products over the dynamics within any invariant subspace, and moreover they do not possess 
normal parameters 

The authors |^,^ have recently shown that dropping these assumptions can lead to the presence of a number of 
novel types of dynamical behaviour, including a new type of intermittency, referred to as in-out intermittency. The 
presence of this type of intermittency has also been found in different distinct nonlinear dynamical systems [^,^. 
Furthermore, there have been interesting developments concerning the study of other phenomena - e.g. riddling ~ in 
these more general settings ^. 

To characterise in~out intermittency, it is best to contrast it with on-off intermittency, as they both can occur in 
systems with invariant submanifolds. To begin with, it is useful to bear in mind that even though on-off intermittency 
can occur in non-skew product settings, all its necessary ingredients can be satisfied in skew product settings. In-out 
intermittency, on the other hand, requires the absence of skew product structure for its existence. 

Briefly, we say that an attractor A exhibits in-out intermittency to the invariant submanifold Mi, if the following 
are true Q| : 

1. The intersection Aq = A Ci Mi is not necessarily a minimal attractor, i.e. there can be proper subsets of Aq 
that are attractors (for on-off intermittency Aq is assumed to be minimal). This means that there can be 
different invariant sets in Aq associated with attraction and repulsion transverse to Aq, hence the name in-out. 
These growing and decaying phases come about through different mechanisms within Mi. If the system has 
a skew-product structure, in-out intermittency reduces to on-off intermittency jij. Fig. |l| shows a schematic 
representation of a typical trajectory for an in-out process near M/. 
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FIG. 1. Typical trajectory of an in-out intermittent solution close to the invariant submanifold Mj, with the two compo- 
nents, the "in" phase and the "out" phase. In the invariant submanifold Mi we may have two or more invariant sets, one of 
which is transversely stable and chaotic but non-attracting in Mi and another which is transversely unstable and is a periodic 
attractor in Mi. The injection mechanism, the "in" phase, is quite irregular and can be modelled by a random walk towards 
Mi, while the expelling mechanism, the "out" phase, can be modelled by a growing exponential spiral away from A//. Note 
that the invariant sets in Mi are represented as points only for clarity. 

2. The minimal attractors in the invariant submanifold are not necessarily chaotic (as for on-off intermittency); they 
are very frequently periodic or equilibria. Furthermore, the trajectory remains close to one of these attractors 
during the moving away or "out" phases, with the important consequence that during these "out" phases the 
trajectory can shadow a periodic orbit, for example, while drifting away from Mi at an exponential rate [Q (see 
also §). 

3. The asymptotic scaling of the probability distribution of the duration of laminar phases in the in-out case can 
have two contributions: 

in— out 

Pn ~ an-'/^e("^"^, +7e^-^"^' ^ h + h, (1) 

on— off 

where a>0,(5>/3>0 and 7 > are positive real constants depending on the bias of the random walk 
modelling the "in" phase and the probability of leaking into the deterministic "out" phase (see for details). 
The term Ii corresponds to biased on-off intermittency, while the extra term I2 can cause an identifiable shoulder 
to develop at large laminar sizes n which can help to statistically distinguish in-out from on-off intermittency. 

The authors in Q were motivated by a numerical exploration of a two dimensional map and explored the statistics 
by means of a Markov chain model. Our aims in this paper are twofold. Firstly, we demonstrate the occurrence of 
in-out intermittency in dynamical systems generated by ordinary differential equations (ODE) as well as by partial 
differential equations (PDE). The latter are especially of interest, since they are in principle infinite dimensional and 
also because few examples of intermittent behaviour and their scalings have been shown concretely to occur in such 
models (see e.g. ^). Secondly, by choosing as our models the mean-field dynamo models the occurrence of this 
type of intermittency could be of interest in understanding certain features of solar and stellar variability, and in 
particular we expect that due to its generic features, it may well appear in more detailed and accurate models of solar 
and stellar variability. 
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II. IN-OUT INTERMITTENCY IN MEAN-FIELD DYNAMO MODELS 



Mean-field dynamo models have been employ ed extensively in order to study various aspects of the dynamics of 
solar, stellar and galactic dynamos (e.g. Their rather idealised nature has been criticized by a number of 



authors (see e.g. |12 ). However, such models are thought to capture some of the essential physics of the turbulent 
processes and reproduce many important dynamical and statistical features of the full three dimensional magneto- 
hydrodynamical models (see e.g. and also 

The standard mean-field dynamo equation is given by 

— = V X (u X B + aB - ?7tV X B) , (2) 

where B and u are the mean magnetic field and mean velocity respectively and the turbulent magnetic diffusivity rjt 
and the coefficient a arise from the correlation of small scale turbulent velocities and magnetic fields ||^ . 

In axisymmetric geometry, equation (|^) is solved by splitting the magnetic field into poloidal and toroidal compo- 
nents, B = Bp -|- B^, and expressing these components in terms of scalar field functions 

Bp = V X A{r, e, t)4>, B^ = B(r, 9, t)$, 

in spherical polar coordinates {r,6,<l)). Equation (|^) can then be expressed in terms of equations for the scalars A 
and B, 
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where V ■ A — and we consider a purely rotational velocity u = uigr^ sm9(l). Nondimensionalisation of these 
equations in terms of a length R and a time /rjt produces the convective and rotational magnetic Reynolds numbers 
Ca = oioRMt and C^^ = UoR^/rjt, where ao and ljq are typical values of a and |a;|. 

Solutions to these equations are often considered in the atu limit where the terms in a can be ignored in the equation 
for B, giving a single dynamo parameter D = CaC^ on rescaling. This reflects the fact that, in stellar convective 
zones, rotational shear produces toroidal flux much more effectively than the processes represented by the a terms, 
whilst in the full equations (the so called a^u limit) we retain both Ca and Cuj as two control parameters. 

The equation (||) gives a kinematic dynamo, since the velocity field u is prescribed. As this equation stands there 
is no mechanism to limit the growth of the magnetic field a nonlinear saturation mechanism is usually supplied by 
making a depend on B. This can be done by supplying a closed functional form representing a fixed approximation 
of the nonlinear effect (c.f. p^ , p^ ), or more dynamically, by supplying an auxiliary equation for a (c.f. [ p^t and 
references therein). 

In the following we consider two cases arising from two separate studies h the above PDE model in the a^i 



3 captions) as well as a finite order 



limit with two different algebraic forms for q;(B) (c.f. (iJJT^ and Figs. || and 
truncation of it in the aco limit but with a time dependent form of the a effect in one spatial dimension (this can be 
obtained by averaging (||) over r) and using a spectral expansion . This ODE model possesses a second (alongside 
D) control parameter, the magnetic Prandtl number — vt/vt: where i^t is the turbulent kinematic viscosity, which 
arises from the time dependent equation for a. This model is given by 

^ = -^^A, + f(B._i+B,+i) (4) 

N N 

i=i fc=i 
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N 

N N 
j=l k=l 

where Ai, Bi and Ci are spatially independent coefficients of the spectral expansions of the scalar fields A^ B and a 
respectively, F, H and G are coefficients expressible in terms of i,j and fc, N is the truncation order and D and v are 
the parameters defined above. The detailed derivation of these equations together with a phenomenological study of 
their dynamics is given in ||^. 

We note that the main ingredients necessary for the occurrence of in-out intermittency are present in both these 
models. Both are axisymmetric and possess invariant submanifolds. More precisely, the truncated model (^) with N—A 
is a 12-dimensional system of ODEs with two 6-dimensional symmetric and antisymmetric invariant submanifolds 
given by 

Ms = {0, Si, 0, yl2, 0, C2, 0, S3, 0, ^4, 0, C4}, (5) 



dt 
dt 



Ma = {Ai, 0, 0, 0, B2, C2, A3, 0, 0, 0, S4, C4}, (6) 

respectively. Similarly the PDE model (|^) possesses two invariant submanifolds, the antisymmetric and symmetric 
invariant submanifolds which are given by 

Ma : A{9)^A{-9), B{e) = ~B{~9), (7) 



Ms : A{e)^A{-0), B{e)^B{~ 



(8) 



respectively, where 9 is the latitude. 

If one separates the poloidal and toroidal scalar field components into symmetric and antisymmetric parts then the 
dynamic evolution for the symmetric (antisymmetric) components has contributions from antisymmetric (symmetric) 
counterparts. This means that these equations are of non-skew product type. For the ODE system (Q) this can be 
readily seen by noting that the evolution equation for each component in Mg {Ma) contains components from Ma 
(Ms). For the PDE models, we first note that u in equation (||) is not dynamical: it is prescribed and therefore can 
be viewed as a part of the initial conditions. The non-skew product nature of the PDE models follows in a similar 
way to the ODE models, bearing in mind the form of the equation (^ and those of the invariant submanifolds (|^) 
and d). 

Finally the control parameters D and v appearing in the ODE model (^) are generically non-normal as they enter 
the equations for Ai and Ci for all i. Similarly this is also true for the control parameters Ca and Ci^ in the case of 
the nondimensionalised version of the PDE equations (^). 

In this way, both models possess all the necessary ingredients for the occurrence of in-out intermittency. 
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FIG. 2. First and second Lyapunov exponents for the attractor of the ODE model restricted to Ma with N=4 and 
v — 0.47, clearly showing the presence of wtndows of periodicity. Note that there is always one null Lyapunov exponent due to 
the system of equations (^) being autonomous. 

Given that the ODE models are more transparent, we first demonstrate the presence of in-out intermittency in 
the truncated system (0) with A^=4. For this model, in~out intermittency occurs for parameter values for which the 
system of equations (^restricted to Ma is within a window of periodicity (c.f . ||l8| ) . Fig. |2| depicts the presence of 
such windows for the ODE model (^. The presence of such windows is supported by a conjecture of Barreto et al. 
psf , according to which for chaotic systems with k positive Lyapunov exponents and m control parameters, with 
m > k, there is a dense set of nearby parameter values at which the attractors are periodic. This implies that for 
our system (^, for each parameter value at which there is a chaotic attractor in Ma there are parameter windows 
arbitrarily close for which the attractor is periodic. 

Fig. 1^ shows an example of in-out intermittency in this system at parameter values D ~ 177.700196 and i> — 0.47. 
We note that even though the interval reported here over which in-out occurs is small, nevertheless there are likely 
to be other intervals (according to the conjecture of Barreto et al., possibly an infinite number of them) over which 
this happens. 

The top panel shows the periodic orbit in the antisymmetric invariant submanifold. Ma, which the projection of the 
trajectory of the full system shadows clearly (second panel). This shadowing or intermittent periodic locking of the 
tangential variables occurs within the laminar phases (third panel) where there is a simultaneous exponential growth 
(hence the name "out" phase) of the amplitudes of the transverse variables through several orders of magnitude 
(bottom panel). This last panel also shows the "in" phases, which can be modelled as a biased random walk taking 
the trajectory into the invariant submanifold. 
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FIG. 3. In-out intermittency in the ODE model (^ with N='i and parameter values D — 
energy and parity are given by _B = E"^ + E^ and P = {E^ — E^)/E respectively, where E^ 
and symmetric parts of the magnetic field energy with respect to the equator ("antisymmetric" {P = —1) and "symmetric" 
[P — +1)). The top panel shows the evolution of an initial condition in Ma and the other panels a nearby initial condition 
not in Ma- In these panels, we have taken a Poincare section at Ai—O for clarity and comparison. 
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FIG- 

model (^) for the case considered in Fig. using a time step r = 10~^n. The shoulder at large laminar phases (which identifies 
the influence of I2 and is a characteristic of in-out intermittency) is easily discerned. 



To substantiate this further, we also calculated the scaling of the probability distribution of the duration of laminar 
phases and this is shown in Fig. ^ This is compatible with the predicted scaling (|^), possessing both a n~'^/^ section, 
at small laminar phase sizes, as well as a noticeable shoulder at higher laminar phase sizes, the latter being a distinctive 
signature of the in-out intermittency. 
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FIG. 5. In-out intermittency in the axisymmetric PDE mean-fleld dynamo model (|3|). The parameters used were ro = 0.4, 
Ca ~ 1.942, — —10^, / = 0.0, with the usual algebraic form of a = ao/(l + B^) (see for details of the parameters). To 
visually enhance the periodic locking we have time sampled the series in the two upper panels. 
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FIG. 6. In-out intermittency in the axisymmetric PDE mean-field dynamo model (|3|J 
Ca ~ 1-5, Cui = —10"", / = 0.7, together with an algebraic form of a due to Kitchatinov 
as in Fig. n. 



The parameters used were ro — 0.4, 
15|. The two upper panels are shown 



These signatures, namely the periodicity of the attractor of the system restricted to the invariant submanifold, the 
periodic locking and the exponential growth of the "out" phases together with the compatibility with the scaling (|l|) 
clearly show the occurrence of in-out intermittency in the truncated ODE dynamo systems. 

To demonstrate the occurrence of in-out intermittency in the PDE case (which as shown above also possesses 
all of the required ingredients), we integrated equation in parameter regions suggested by fl^ ], using the code 
described in and implemented by Figs. ^ and ^ give examples of in-out intermittency in these PDE models 
pof . As can be seen, this behaviour can occur with the invariant submanifold being either antisymmetric (Fig. |^) 
or symmetric (Fig. |^). Again, in addition to the presence of periodic behaviour in the system restricted to the 
invariant submanifold (top panels), these figures clearly show the presence of locking during the "out" phases (second 
panels) with an exponential growth of the energy of the transverse modes through several orders of magnitude (bottom 
panels). This behaviour mirrors very closely the truncated ODE model shown in Fig. ^ as well as that expected to 
occur from the theory Q. To substantiate this further, we again looked at the compatibility of the scaling for the 
distribution of the laminar phases with the theoretical scaling given by (|^) . Despite the greatly enhanced numerical 
cost of integrating the PDE equations long enough to obtain convergence to the scaling law, we have been able to 
establish agreement in this case as can be seen in Fig. ^. Together, these signatures clearly demonstrate the occurrence 
of in-out intermittency in these PDE dynamo models. 
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FIG. 7. Scaling of the probability distribution of the duration of laminar phases for the axisymmetric PDE mean-field 
dynamo model (^) for the case considered in Fig. ^, using a time step r = 10~*n. The shoulder at large laminar phases (which 
identifies the infiuence of I2 and is characteristic of in-out intermittency) is easily discerned. 



III. DISCUSSION 



By establishing the main ingredients necessary for the occurrence of in-out intermittency as well as checking 
the predicted corresponding phase space signatures and predicted scalings, we have concretely demonstrated the 
occurrence of this type of intermittency in both ODE and PDE models. This type of intermittency requires for its 
existence the non-skew product feature, the generality of which makes the occurrence of this type of intermittency of 
potential interest. 

The models chosen here are mean-field dynamo models, which despite their approximate nature are thought to 
capture many features of magnetic activity in solar-type stars. An important observed feature of variabilities in 
solar-type stars is the presence of dynamical behaviour with different statistics over different time intervals due to the 
occurrence of the so called grand minima during which the amplitude of the magnetic activity is greatly diminished. A 
number of scenarios have been suggested in order to explain these phenomena (see e.g. [p2|). Within the deterministic 
framework, intermittency (and multiple intermittency p^ ) has been put forward as a possible mechanism. A 
number of studies have found intermittent types of behaviour in such models (e.g. and references therein). The 
concrete demonstration of in-out as well as other forms of intermittency are of potential importance in this regard as 
they demonstrate the possible types of dynamical variability that can occur in such settings. 

We thank Axel Brandenburg and David Moss for helpful conversations. EC is supported by a PPARC postdoctoral 
fellowship, PA was partially supported by EPSRC grant GR/K77365 and RT benefited from PPARC UK Grant No. 
L39094. 
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